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We characterize the operation of semiconductor micro-ring lasers in an excitable regime. Our 
experiments reveal a statistical distribution of the characteristics of noise-triggered optical pulses 
that is not observed in other excitable systems. In particular, an inverse correlation exists between 
the pulse amplitude and duration. Numerical simulations and an interpretation in an asymptotic 
phase space confirm and explain these experimentally observed pulse characteristics. 



PACS numbers: 42.65.Pc, 42.55. Px,42.60.Mi 

I. INTRODUCTION 

Excitability — the ability of nonlinear systems to 
fire large well-defined output pulses when a threshold is 
crossed — has been investigated in a broad range of scien- 
tific areas, including physics, chemistry, biology and neu- 
roscience p~HI] • The strong interdisciplinarity of excitable 
systems has attracted the attention of theoreticians in- 
terested in their universal properties. These theoretical 
studies disclosed that excitability takes place when a sep- 
aratrix is crossed in the phase space of the system, and 
led to a comprehensive classification of excitable systems 
based on bifurcations taking place in their phase space 
[31 |B]. Considerable theoretical insight was further gained 
into the role played by networks of excitable units in 
neuroscience [8]. Moreover, the investigation of excitable 
units is of high practical relevance as networks of spiking 
neurons have proved to be computationally superior to 
other neural networks @]. 

In the field of photonics, excitability is widely stud- 
ied both theoretically and experimentally [9HT7] and a 
broad search was started for optical excitable units that 
could be deployed in optical neural networks. In the last 
decade, lasers with saturable absorber [9l[T0], optically 
injected lasers [HI US], lasers with optical feedback [15] 
or VCSELs with opto-electronic feedback [T7] have all 
been proposed as optical excitable units. 

Most often excitable behavior in optical systems has 
been shown to occur close to a fold bifurcation and a 
homoclinic bifurcation of a stable limit cycle [5rfl2l HB] . 
In these systems, the response to noise is governed by 
the presence of an accessible saddle point S embedded 
in the separatrix. Pulses are activated by noise-induced 
fluctuations which connect the resting state to S and 
the large deterministic excursion takes place along the 
unstable manifold of S. In this scenario, excitability is 
therefore possible in the limit of vanishing noise intensity. 

In a recent paper [IS], we proposed that excitabil- 
ity can take place in systems with weakly broken Z2- 
symmetry, which includes optical units such as semicon- 
ductor micro-ring lasers [19j [20] and micro-disk lasers 
[2"T] . It can be shown that in this class of systems the ac- 



cessible saddle is not embedded in the separatrix for ex- 
citability. Therefore, the unstable manifold of S does not 
participate in excitability and the generation of pulses 
cannot be initiated in the vanishing noise limit. Clear 
evidence of such excitability was shown in semiconductor 
ring lasers with weakly broken Z2-symmetry in the pres- 
ence of a non-vanishing, finite intensity of spontaneous 
emission noise [18 . 

In this paper, we address that micro-ring lasers 
(MRLs) in the excitable regime show a notable degree 
of variation in the amplitude and width of the excited 
pulses. We characterize both theoretically and experi- 
mentally the particular pulse properties when the trig- 
gering occurs through optical noise in the system. The 
main result is that a clear inverse correlation is observed 
between the amplitude and width of the excited pulses. 
In Section [n] we experimentally investigate the excitable 
behavior of a MRL in the time domain. A characteri- 
zation of the stochastic properties of excitability is car- 
ried out to reveal an inter-spike-interval (ISI) distribution 
which is exponential for large waiting times, but diverges 
from the Kramers form for intervals below approx. 50 ns. 
In Sections |III| and |IV[ we use a general rate-equation 
model and an asymptotic model for excitability in MRLs 
as introduced in Ref. [HI [22] to explain the experimen- 
tally observed features. Finally, in Section [Vj we discuss 
the main results presented and their generality. 



II. EXPERIMENTS 

A. Device and Setup 

The experiments have been performed on an InP-based 
multi-quantum- well MRL with a racetrack geometry [see 
Fig. [it a)]. The optical power is coupled out of the ring 
cavity by directional coupling to bus waveguides which 
are integrated on the same optical chip. The use of two 
bus waveguides in the micro-ring design allows for four 
independent input/output ports A-D which can be ac- 
cessed with optical fibers. The device is mounted on a 
copper mount and is thermally controlled by a Peltier 



2 



(a) 

I : Waveguide current 
I : Ring current 



0.05 



Lensed 
Fiber 




(b) 



SOA 






SRL 








read-out port A 



FIG. 1. Experimental MRL set-up. (a) depicts the actual 
MRL device with its contacting. The ring itself is pumped 
with a current I p . Four waveguide contacts are depicted of 
which only port A is biased with I w . (b) shows the overall 
setup where the read-out and the SOA are separated from 
each other by a circulator C. 
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element which is stabilized with an accuracy of 0.01°C. 

The MRL chip is deployed in a setup as shown in Fig.[l] 
In our experiments only port A of the four ports is used 
for collecting output power from the MRL as well as in- 
jecting optical noise in the micro-ring cavity. We use a 
directly biased semiconductor optical amplifier [SOA in 
Fig. [ljb)] operating in the C-band to generate a con- 
trollable amount of noise through amplified spontaneous 
emission which is injected in the ring through a lensed 
optical fiber coupled to port A of the MRL. The use of a 
circulator [C in Fig. [l](b)] allows us to read output from 
port A with the same lensed fiber. 

Electrical contacts have been applied to the bus waveg- 
uides that can be independently pumped. They allow to 
amplify on-chip the signal emitted by the ring. More in- 
terestingly, the presence of a contact allows us to contin- 
uously break the symmetry of the device in a controlled 
way by making both the strength and the phase of the lin- 
ear coupling between the CW and the CCW mode asym- 
metric. Using the fiber's facet as a mirror, we are able to 
reflect power from one mode (e.g. CCW) back into the 
waveguide and finally to the counter-propagating mode 
in the ring. The application of a direct bias current I w 
on the waveguide's electrode has two main effects. First, 
it controls the power that is coupled to the CW mode; 
second, it affects the optical length of the waveguide via 
carrier-induced refractive index changes. Therefore, the 
phase of the reflected signal can be controlled by tun- 
ing the bias current on the waveguide and changing the 
position of the fiber facet which is piezo-controlled. We 
have previously demonstrated in Refs. [18, 22 that such 
a breaking of the Z2-symmetry of the micro-ring leads 
to excitable behavior when the optical noise intensity is 
large enough to make one of the states metastable. We 
analyze the output power of the CCW mode (port A) 
using a fast photodiode (with a bandwidth of 2.4 GHz) 



FIG. 2. (a) Experimental time series of port A demon- 
strating the excitable behavior of the MRL operating at 
I p — 47.45mA. The current on the waveguide of port A 
is I w = 13mA and the current in the external SOA is 
Isoa = 700mA. (b) and (c) show examples of a single and a 
double excited pulse taken from the same time series. 



connected to an oscilloscope (with a maximum sampling 
rate of 20 Gs/s). 



B. Excitability 

In the temperature range of operation, the trans- 
parency current density of our semiconductor material is 
~ 1.0kA/cm 2 , which leads to a transparency current of 
~ 24mA for the MRL and ~ 4mA for the waveguide. The 
MRL device reaches threshold at 35mA and excitability 
is observed between 42mA and 48mA. A typical time 
series revealing excitable behavior of the micro-ring is 
shown in Fig. [2ja) . The device operates most of the time 
in the CW unidirectional mode. However, pulses in the 
CCW direction can be regularly observed. An example of 
an excited pulse with a pulse duration of ~ 7ns is shown 
in Fig. Kb). 

In contrast to what is expected for a typical excitable 
system, not all excited pulses have the same amplitude 
and duration, see e.g. Fig. [2ja). We clearly observe a 
degree of variation in the amplitude of the pulses (which 
is not a consequence of undersampling since it persists for 
higher sampling rates). In the same way, a distribution 
of pulse-durations is observed in the experimental time 
series and the pulse-durations seem to be spread between 
a minimum of 7ns and a maximum of 20ns. Such a spread 
is significantly larger than the sampling time. 
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FIG. 3. Inter-Spike-Interval (ISI) distribution for pulses emit- 
ted by the symmetry broken MRL for device parameters as 
given in Fig. [2] The logarithm of the number of events N in 
each ISI bin is shown, (a) The distribution of 'long' ISIs is 
well- fitted by an exponential curve (see gray solid line). The 
gradient 1/Ti = 1.6/xs -1 indicates a characteristic time-scale 
of ~ 0.63 /us. (b) The distribution of 'short' ISIs reveals a 
different time constant Tb ~ 18.4ns. The inset in panel (a) 
shows the dependence of the 'long' ISI time-scale Ti (with er- 
ror bars given by the standard deviation from the mean value 
Ti) on the current on the external SOA (Isoa) which is used 
as noise source. 



C. Inter-spike-interval distribution 

In order to quantitatively describe the ISI, we have 
constructed the ISI-distribution from the experimental 
time-series. A threshold is put at 60 % of the maximum 
pulse height recorded in the time series. Only pulses 
above this threshold are taken into account in the cal- 
culation of the ISI distribution. We have verified that 
calculating the ISI distribution for other values of the 
threshold gives the same qualitative results as long as it 
is not chosen too low as one will then start to measure 
noisy excursions instead of actual excited pulses. An ex- 
ample of such an ISI-distribution is shown in Fig.^a) for 
the MRL operating at a current of 47.45 mA. The tail of 
the distribution is well fitted by a typical Kramers' expo- 
nential distribution P (r) oc exp [— t/7\], which indicates 
that the excited pulses are activated by noise-induced 
crossing of a potential barrier. The time constant Ti ob- 
tained from the fitting of Fig. [3|a) is Ti = 0.63 fjs. 

A significant deviation from the Kramers' law is, how- 
ever, observed for ISIs shorter than 50 ns, as shown in 
Fig. [3^b) . The distribution of such 'short' ISI is strongly 
peaked around approx. 20ns and decreases abruptly with 
a characteristic time T 2 = 18.4ns. A similar devi- 
ation from the Kramers' law was previously reported 
in the residence-time distribution of MRLs [53] and in 
other optical excitable systems such as lasers with optical 
feedback close to the low frequency fluctuations regime 

El |25] . 

In order to confirm and quantify the noise-activated 
origin of excitability, we have investigated the depen- 
dence of the average Inter-Spike-Interval T\ as function 
of the bias current on the ring I p and the external noise. 
The external noise source is the external semiconductor 
optical amplifier [see Fig. [ljb)] and it is biased at dif- 
ferent DC currents. We assume a linear relationship be- 



FIG. 4. Distribution of (a) the pulse amplitude and (b) the 
pulse width for I p = 47.45mA, Isoa = 700mA, I w = 13.0 
mA. The number of events N in each bin is plotted. The 
gray solid line shows the corresponding Gaussian distribution 
as predicted by the theory presented in Section [TV] 

tween the SOA current and the noise intensity coupled 
into the MRL. In the inset of Fig.|3ja) the dependence of 
Ti versus the SOA current is shown. The linear relation 
between log(Ti) and the inverse of the SOA current is 
consistent with a Kramers rate across a potential barrier 
[26] . We have also verified that Ti increases with the 
bias current I p on the ring, which is consistent with the 
theoretical result that the activation energy of the MRL 
increases with the bias current [27] . 

In contrast, we have checked that the 'short' ISI has a 
much smaller dependence on the noise strength and re- 
mains approximately constant when changing the Isoa 
and I p . The origin of the non-Arrhenius distribution for 
these 'short' ISIs is instead related to deterministic dy- 
namics. This will be discussed in more detail in Sections 

nnna 



D. Amplitude-width distribution 

In order to further characterize the properties of the 
excited pulses, we build the distribution of the pulse du- 
ration, defined as the full- width-half-maximum (FWHM) 
and the distribution of pulse amplitudes. Typical exam- 
ples of the amplitude and width distributions are shown 
in Fig. H] 

A distribution of pulse amplitudes is evident in 
Fig. |4]ja), which is consistent with the pulse amplitude 
modulation present in Fig. [2] Such distribution of ampli- 
tudes was not observed in other optical excitable systems 
(see for instance [HI HH [2S] ) an< ^ * s no * consistent with a 
standard excitability scenario where the system performs 
a large phase space excursion following the unstable ma- 
nifold of a saddle structure [5J[7] for which a much sharper 
distribution is expected. The width of the excited pulses 
is distributed asymmetrically around an average value of 
18 ns and reveals a non-zero probability for pulses as long 
as 40 ns. 

In order to shed light onto the dynamical origin of such 
distributions, we have measured the correlation between 
pulse amplitude and width for different values of the bias 
current on the ring. We minimize the noise-contributions 
to the pulse amplitude by performing this measurement 
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FIG. 5. Measured width and amplitude of the excited pulses 
for two different values of the bias current I p on the ring. 
Isoa = 0mA and I w = 13.0 mA. (a) I p = 42.82mA and (b) 
I p — 44.49mA. A linear fit of the data points is shown by the 
gray, solid line. The correlation coefficient C is approx. 0.7 in 
both cases. 



without the external SOA. The results are shown in Fig. [5] 
and quantified using a correlation coefficient C: 
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For all used bias currents there is a clear correlation 
between the width and amplitude of the pulses. The 
experiments reveal that pulses with a higher peak power 
are narrower, whereas pulses with a lower peak power 
are wider. For Isoa = 0mA and /„, = 13.0 mA, we 
have measured a correlation C of about 0.7. The linear 
fit between the amplitude and width is indicated by the 
gray, solid line. We remark that this correlation shown in 
Fig. [5] has been obtained by defining the pulse-width by 
its FWHM; we have verified that other ways to quantify 
the pulse duration (not shown) lead to consistent results. 

To understand the origin of this amplitude-width cor- 
relation of the excited pulses, we will numerically analyze 
the pulse properties in Section |IV| To allow for a clear 
presentation of the results, in the next section, we shortly 
retake our approach towards modeling the time evolu- 
tion of a single-longitudinal, single-transverse mode in a 
micro-ring cavity using a general rate-equation model for 
MRLs and an asymptotically reduced model. In particu- 
lar, we explain the excitability scenario for MRLs in the 
asymptotic phase-space. For more information on the 
rate-equation model, we refer to Refs. [2"2l HZ51 HZtj] . 



III. MODELING AND ORIGIN OF 
EXCITABILITY IN SEMICONDUCTOR RING 
LASERS 

In order to model the MRL operating in single- 
longitudinal and single-transverse mode, we use a rate- 
equation model for the the evolution of the slowly vary- 
ing amplitudes of the counter-propagating modes E cw ccw 
and the carrier inversion N 1221: 



Symbol Physical meaning 


Simulation value 
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Field decay rate 


100ns" 1 


7 


Carrier inversion decay rate 


0.2ns" 1 
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Linewidth enhancement factor 


3.5 


/' 


Renormalized bias current 


1.65 
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Self-saturation coefficient 


0.005 
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Cross-saturation coefficient 


0.01 


k 


Coupling amplitude 


0.44ns" 1 


Ak 


Coupling amplitude asymmetry 0.044ns 1 


<t>k 


Coupling phase 


1.5 




Coupling phase asymmetry 






TABLE I. Summary of the physical meaning of the parame- 
ters in the rate equations |3]l-([5| and their typical values used 
throughout this article, unless stated otherwise. 



df 



dE r 



dt 

dN 
~dt 



k(1 + ia) [g cw N - 1] E cw 

-(k - Ak/2)e^~ A ^^E ccw + £ cw , (3) 

k(1 + ia) [g ccw N - 1] E ccw 

-(k + Ak/2)e^ +A ^E cw + U w , (4) 

7 [A* — N — g cw N\E cw \ 2 - g ccw N\E ccw \ 2 } (5) 



where g cw = I- s\E cw \ 2 - c\E ccw \ 2 , g ccw = l-s\E ccw \ 2 - 
c\E cw \ 2 is a differential gain function which includes phe- 
nomenological self (s) and cross (c) saturation terms, k 
is the field decay rate, 7 is the inversion decay rate, a is 
the linewidth enhancement factor of the semiconductor 
material and /i is the renormalized injection current with 
fj, w at transparency and fi w 1 at lasing threshold. 
For the device such as the one discussed Sec. |H} \i = 
corresponds to 25mA and /1 = 1 to 35mA, yielding the 
approximate relation I p w (10/i + 25)mA. 

Asymmetric linear coupling terms are present in 
Eqs. Q-Q between E cw and E ccw , which model a 
backscattering of power from one mode to the other. 
Such intrinsic backscattering originates from reflections 
at the directional coupler or at the chip facets, and is in 
general asymmetric due to unavoidable imperfections in- 
troduced during device fabrication. Moreover, asymme- 
tries in backscattering introduced externally in our set-up 
(see Fig.[T]), such as reflections at the fiber tip at one side 
of the chip, are also lumped into these coupling terms. 
The parameters k and <^ represent respectively the aver- 
age coupling amplitude and phase, whereas the coupling 
asymmetry is described by the symmetry breaking terms 
Afc and Aipk- Unless mentioned otherwise, throughout 
this manuscript we will use the parameters shown in Ta- 
ble m 

Noise terms are introduced in Eqs. ([3])-([5]) as complex, 
Gaussian, zero-mean stochastic terms £ C w,ccw described 
by the correlation terms + r)£?(t)) = 2DN8 1] 5(t) 1 
where i,j ={cw, ecw} and D is the noise intensity |30j . 
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FIG. 6. The phase space topology of an asymmetric MRL. 
CW, CCW and S are respectively two quasi-unidirectional 
states and a saddle state of which the stable and unstable 
manifolds are displayed. The basin of attraction of the CW 
(CCW) quasi-unidirectional state is colored white (gray). Pa- 
rameters are chosen as in Table [U 



All simulations are performed in this rate equation 
model ([3])-([5]). However, we will interpret our results in a 
two-dimensional phase-space corresponding to a reduced 
MRL model that has been proven to be a powerful tool to 
model the slow time dynamics in MRLs [15112211251151] . 
This two-dimensional model is valid on time scales slower 
than those of the relaxation oscillations. The two vari- 
ables 9 € [— 7r/2,7r/2] and ip £ [0, 2tt] are defined by 



9 = 2 arctan 



E, 



2' 



(6a) 
(6b) 



9 is a measure for the power partitioning between the 
counter-propagating modes and ip is the phase difference 
between the corresponding electric fields. 

To describe the excitability scenario in MRLs, we con- 
sider a ring laser operating unidirectionally whose sym- 
metry has been slightly broken. The corresponding phase 
space is shown in Fig. [6j Two counter-propagating unidi- 
rectional stable attractors are present in the MRL — the 
CW and the CCW mode — and are depicted here in the 
two-dimensional (#, ip) phase space. The white and gray 
regions indicate the basins of attraction of the CW and 
the CCW mode. They are separated by the stable man- 
ifolds of a saddle state indicated by S. 

The presence of an asymmetry parameter manifests 
itself in the different stability of the modes, leading to 
different sizes of their basins of attraction. In particular, 
the distance between two branches of the stable manifold 
of S is affected by the symmetry breaking and can be 
made arbitrarily small by controlling the parameters of 
the system [22]. When noise is present in the system, a 
diffusion length-scale Ld appears, which depends on the 
noise intensity D. The onset of excitability in MRLs is 
regulated by the interplay between Ljj and the distance 
between the folds of the stable manifold of S. 

Assume that the MRL is operating in a regime such 
that Ljj is small compared to the size of the basin of 
attraction of CW and large compared to the distance 



between the branches of the stable manifold of S. The 
MRL will spend most of the time in the vicinity of the 
CW stable state; however a rare large fluctuation may 
move the system to the boundary of the basin of attrac- 
tion of CW. When this happens, the system will cross this 
boundary with an overwhelming probability by crossing 
both branches of the stable manifold of S and thereafter 
perform a large deterministic excursion leading to the 
emission of a CCW pulse. 

This excitability scenario is different from the one most 
frequently encountered in optics, where pulses are initi- 
ated stochastically by crossing only one branch of the 
stable manifold of a saddle and completed largely deter- 
ministically by following a branch of the unstable mani- 
fold of the same saddle back to the initial quiescent state 
[HHHJ HH US] • Even in the rare situations where the sys- 
tem was shown to be both bistable and excitable at the 
same time, the excitable excursions are still completed 
by following the unstable manifold of the saddle point 
[32] . In MRLs, however, this scenario is forbidden by 
the residual Z2-symmetry and the unstable manifold of 
S which connects with the metastable CW state [TS] . 

This specific mechanism of excitability is general for 
all systems with weakly broken Z2-symmetry and occurs 
near a homoclinic bifurcation that unfolds from a Takens- 
Bogdanov point. In that homoclinic bifurcation an un- 
stable cycle is created, which later disappears in a fold of 
cycles [22] . Such a sequence of bifurcations leads to the 
folded shape of the stable manifold of the saddle (sec Fig. 
[6]), necessary for the system to be excitable. The unfold- 
ing of the different bifurcations from a Takens-Bogdanov 
point has been characterized in depth, both in systems 
with Z2-symmetry 33, 34 a and in systems where this sym- 
metry is broken [35] . As an example of other systems that 
share the same symmetry as the MRL and which can 
therefore exhibit similar excitable behavior when weakly 
breaking this symmetry, we mention e.g. C02-lasers [36] 
and oscillatory convection in binary fluid mixtures |37j . 
Perhaps the most obvious example with the same cir- 
cular symmetry as MRLs are semiconductor micro-disk 
lasers |21j . In semiconductor micro-disk lasers where the 
two lasing modes are the whispering gallery modes, the 
appropriate rate-equation model is identical to the one 
studied in this work. 



IV. STOCHASTIC ANALYSIS AND 
COMPARISON TO THE EXPERIMENTS 



In this Section, we use direct numerical integration 
of Eqs. ([§-([5) using a stochastic Euler-Heun method 
to characterize the specific features of excited pulses in 
MRLs. More specifically, numerical time-series are col- 
lected and statistics of the ISI, the pulse amplitude and 
the pulse width are built. A projection of the time-series 
on the reduced phase space (9, i/j) is then performed in 
order to validate the topological arguments of Sec. |III| 
explaining the spread in amplitude of the pulses and the 
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FIG. 7. (Color online) Results of simulations of Eqs. q3p- 
®. (a)-(b): D — 1.5 x 10~ 4 ns~\ t obs = 100 ms, (c)-(d): 
D = 2.5 x 10~ 4 ns~\ tabs = 10 ms, with D the noise strength 
and t bs the observation time, (a) and (c): Histogram of the 
trajectories in asymptotically reduced two-dimensional phase 
space. The full (dashed) white lines indicate the stable (un- 
stable) manifolds of the saddle, (b) and (d): Pulse width vs 
pulse amplitude. The white curve indicates the prediction 
from the deterministic reduced model. Parameters are taken 
as in Tabled] 



associated correlated spread in their width, as also ob- 
served experimentally (see Figures 4|5|. 



A. Pulse amplitude and width 

Fig. [7] shows the distribution of the phase space trajec- 
tories corresponding to excited pulses. In order to avoid 
sampling of double pulses, the system is reinitialized to 
the original CW mode when the tail of the pulse sat- 
isfies the condition 9 > 0. Rare excursions that reside 
longer than 20 ns around the metastable CCW state are 
discarded. Figs. [7]ja) and (c) show histograms of the 
collected trajectories projected on the asymptotically re- 
duced two-dimensional phase space (6,1/1), and this for 
two different values of the noise strength D. The full 
(dashed) white lines indicate the stable (unstable) man- 
ifolds of the saddle S. 

It is clear from Fig. (7^a) and (c) that noise-activated 
trajectories cross the excitability threshold at a finite dis- 
tance from the saddle S. Such saddle-avoidance is ex- 
pected in stochastic systems where is not the shortest 
length scale [38] and it is therefore compatible with the 
mechanism described in Sec. IIIII 

After crossing the stable manifolds of S, the trajecto- 
ries spread in the phase space due to diffusion. A dis- 
tribution of pulse amplitudes and pulse widths as ex- 
perimentally observed in Figs. [2] and [4] can therefore be 
expected in MRLs as the deterministic evolution of the 
pulse does not take place along a unique trajectory. In 
the reduced phase space, pulses are a one-parameter fam- 



ily of trajectories that can be parameterized by their ini- 
tial conditions (0i n ,ipin) [see Fig.©. By fixing 9 m = 
and ipin beyond the excitability separatrix, a sampling of 
excited pulses is achieved. The width of each pulse can 
be quantified by the required time to return to the 9 > 
condition; in the same way, the extreme value 9 max can 
be used to quantify the pulse amplitude. 

The correlation between amplitude and width can now 
be understood due to the different velocity-fields at differ- 
ent positions in the reduced phase space. Higher (larger 
\9\) pulses move faster in phase space and are thus nar- 
rower, while lower pulses are slower and consequently also 
wider. The correlation curves extracted in such a way are 
plotted in Fig. ^b) and (d) (see the solid, white line). 
The observed increase in pulse width with the decrease 
in pulse amplitude confirms the experimental trend re- 
ported in Sec. II D i.e. we refer to Fig. [5} A projection of 



the pulses, obtained by numerically solving the stochastic 
full-rate equation model, in the amplitude-width space 
[see point cloud in Fig. W^b) and (d)] provides an extra 
confirmation of this deterministic prediction. 

We observe that the theoretical amplitude- width curve 
is not evenly sampled by the numerical pulses and 
that the numerical data clusters around a middle value. 
For this reason, we argue that the full profile of the 
amplitude-width curve cannot be observed experimen- 
tally in our devices. The position of the point cloud 
in the amplitude-width space is slightly affected by 
noise. The correlation relation holds for the two different 
noise strengths, although it is clear that for lower noise 
strengths [Fig. 7^b)] the pulses tend to be less high in 
amplitude and thus wider than in the case of higher noise 
strengths [Fig. 7fd)]. The fact that the point cloud in 
Fig. gb) and (d) opens up for decreasing pulse ampli- 
tudes is related to the nature of these pulses. These rela- 
tively long pulses tend to wander around the metastable 
CW state or the saddle point. Near the saddle point 
the deterministic trajectory greatly slows down. Hence, 
in this case, the transit times of the pulses are mainly 
determined by the noise, giving rise to a larger spread. 

A more tangible picture is given by considering 
the width distribution and the amplitude distribution 
separately, instead of their mutual projection on the 
two-dimensional (9, ip) phase-plane. These histograms 
are shown in Fig. [8j The pulse amplitude distribu- 
tion can be properly fitted by a Gaussian distribution 
Aexp[—(9 max — B) 2 /C 2 }, indicating that the amplitude 
of the pulse is mainly determined by the magnitude of the 
perturbation and less by the topology of the flow. The 
average pulse amplitude also increases with increasing 
noise intensity, initiating excursions farther away from 
the stable saddle manifold. In contrast, the pulse width 
distribution is asymmetric. Some of the pulses tend to 
erratically wander around the CCW state or the saddle 
giving rise to the tail in the pulse width distribution. 
The average pulse width decreases with increasing noise 
intensity indicating that the flow near the stable mani- 
fold is slower. This also confirms the pulse amplitude- 
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FIG. 8. (Color online) Results of simulations of Eqs. |3|-([5|. 
The histogram of events N/N to t is shown in function of the 
pulse amplitudes (a),(c) and in function of the pulse widths 
(b),(d). D = 1.5 x 1CT 4 for (a)-(b) and D = 2.5 x 10~ 4 for (c)- 
(d). The histogram of the pulse amplitudes can be fitted by 
a Gaussian: y4exp(— {{x — B)/C) 2 ). (a) N to t = 15052 events, 
A = 0.0575, B = 0.9762, C = 0.1237. (c) N tot = 23748 events, 
A = 0.0503,5 = 1.019, C = 0.1508. The vertical gray lines 
show the maxima of the histograms: a width of ~ 9.20 ns is 
found in (b) and ~ 7.99 ns in (d). Parameters are taken as in 
Table. U 



width trade-off trend of Fig. [7] Similar experimental 
histograms have been shown in Fig. |4j Fig. Qa) showed 
that the spread in pulse amplitudes is consistent with a 
Gaussian distribution, while the spread in the width of 
the pulses in [see Fig.Qb)] is distributed asymmetrically. 

Finally, in Fig. [9] we show typical pulse shapes as the 
pulse width varies. Figs. |9][a)-(e) show a random col- 
lection of pulses with widths in the following intervals: 
5-6 ns, 8-9 ns, 11-12 ns, 14-15 ns, 17-18 ns, respectively. 
The corresponding trajectories in the (0,ip) phase space 
are shown in gray. The faster, narrower pulses start out 
at a relatively large distance from the stable manifold, 
and remain distant from it during the whole pulse tra- 
jectory. This is clearly visible in the pulses in Fig. [9ja)- 
(b). Oppositely, pulses that start out closer to the stable 
manifold slow down and are less high. For this type of 
excursions, if the pulse trajectory comes too close to the 
stable manifold it can get tangled up in the metastable 
CCW state or slowed down near the saddle, explaining 
the formation of possible plateaus in the pulse. In the 
case one gets trapped in the CCW state, the pulse would 
show a plateau at the pulse maximum, while if the pulse 
slows down near the saddle, it would exhibit a plateau 
at the trailing edge of the pulse. In Fig. [9jcl) , many of 
the pulses come very close to the saddle point, result- 
ing in an even more pronounced slow-down, especially at 
the trailing edge of the pulse. Several rare pulse events 
where the system briefly gets stuck in the CCW state are 
depicted in panel (e). Also pulses that do not get close 




FIG. 9. Simulation of Eqs. |3}-|5). Typical pulse shapes cor- 
responding to different pulse widths are shown. Single pulses 
are collected from the time traces in panels (a)-(e) with next 
to each panel (a)-(e) the corresponding trajectory in the (8, if)) 
phase space in gray. Stable (full) and unstable (dashed) man- 
ifolds are drawn in black. From (a) to (e) pulses are shown 
with widths in the following intervals: 5-6 ns, 8-9 ns, 11-12 



ns, 14-15 ns, 17-18 ns. D = 2.5 x 10" 
parameters are taken as in Table. [I] 



ns and the other 



to the CCW state or the saddle point are characterized 
by an asymmetric pulse shape. Such asymmetry due to 
a slowing down of the pulse at the trailing edge (higher 
values of ip) is best visible in panels (b) and (c). 



B. Inter-spike-interval diagram 

We have also determined the ISI distribution by sim- 
ulating Eqs. ([3])- ([5]) during a time t ba = 10ms. The dis- 
tribution of the ISI between all excited pulses is shown 
in Fig. [TUj We have plotted the logarithm of the nor- 
malized number of events log(N/N tot ) in function of the 
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in Figures [7[9j The bin-size has been chosen to be 
ns. It becomes evident from Fig. 



10 'a) that the total 
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FIG. 10. Simulation of Eqs. (|3|-((5| for 10ms showing the 
inter-spike-interval diagram. log2(N/Nt t) is plotted in func- 
tion of the inter-spike-interval. In (a) n — 1.65, D = 
2.5 x 10 -4 ns -1 and the bin-size is chosen to be Ins. In panels 
(b)-(d), a zoom of the peak at shorter ISI is shown with a bin- 
size of 0.05ns. Panels (b) and (c) show the variation of the 
ISI with fj, for D = 2.5 x 10" 4 ns -1 and D = 1.5 X 10" 4 ns" 1 , 
respectively. Panel (d) demonstrates the change in ISI for 
asymmetric noise contributions to the CW (D cm = D) and 
the CCW (D ccw ) mode for fj, — 1.65. The other parameters 
are taken as in Table. [I] 



inter-spike-interval, with N the total amount of events in 
each bin and N to t the total amount of excited pulses. 

In Fig. |lO{a), we have taken [i = 1.65 and D = 2.5 x 
10~ 4 ns -1 corresponding to the parameter set used before 



ISI distribution is the combination of two different time 
scales. The short ISIs and long ISIs are distributed in a 
different way and together represent a strongly Kramers' 
type of behavior. 

The slow time scale can be fitted by an exponential 
curve [oc exp (— t/T\)\, where the average ISI (Ti) is the 
fitting parameter. This slow time scale T\ is typically 
in the order of us and corresponds to the generation of 
the pulses in an Kramers' type noise- activation process 
across the excitability threshold [2B]. The short ISI times, 
however, are a signature of multiple consecutive excited 
pulses due to noise clustering. The possible excitation 
of double pulses such as those experimentally shown in 
Fig. 2|c) depend on the closeness of the stable and un- 
stable manifolds of the saddle S. When these are close 
enough, noise can excite a second pulse before the sys- 
tem can relax to the quiescent state. The presence of such 
excited multi-pulses will show up in the ISI distribution 
as a sharp peak around the average pulse width. The 
numerical results in Fig. |10[ a) confirm the experimental 
observation of both time scales in the experimental ISI 
distribution presented in Section [IT] (see Fig. [3]). 

The presence of these time scales was also found in 
our study of stochastic modc-hopping in the bistable 
regime in Refs. [23J HZ1 ISI]- Both a mode-hop in the 
bistable regime and an excitation beyond the excitabil- 
ity threshold are described by a noise-activated escape, 
corresponding to the slow Arrhenius time-scale. The fast 
non-Arrhenius character of the ISI finds its origin in a 
noise-induced diffusion through both branches of the sta- 
ble manifold, thus initiating another excursion in phase 
space before relaxing to the CW state. Such a noise 
clustering of pulses due to the proximity of the relax- 
ation trajectory of the excited pulse and the excitability 
threshold has been observed in several other noise driven 
excitable systems, such as e.g. lasers with optical feed- 
back [211I2S], quantum-dot lasers with optical injection 
[39l 140] and neurons of the Hodgkin-Huxley type [41] . 

In panels (b)-(d), we study in more detail the ISI dis- 
tribution of the faster time scale for different values of the 
pump current /i, the noise strength D and for asymmet- 
ric contributions of the noise to both counter-propagating 
modes. The bin-size is taken to be 0.05ns. 
Figs. [Toj^b)-(c) show the ISI distribution for varying val- 
ues of the current ji at a fixed noise strength D = 
2.5 x 10~ 4 ns -1 and D = 1.5 x 10~ 4 ns -1 , respectively. 
One can notice that in all cases the maximum amount of 
events are located around 5ns, which corresponds roughly 
to half of an excitable excursion [32]. For slightly higher 
ISI (5 — 10ns) a dip in the ISI distribution is observed. 
Such a dip is typical for noise clustering [2H HI] and is 
due to the nature of the relaxation trajectory of the ex- 
cited pulse. In particular, it becomes evident from Fig. 
[lOjb)-(c) that the dip becomes less pronounced for de- 
creasing values of the current fi due to the fact that the 
relaxation to the stable node occurs increasingly slowly 
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for lower values of fj,. One can also notice that the entire 
ISI curve moves up for decreasing /i, which is a logical 
consequence of the decreasing depth of the potential well. 
It is interesting to note that the basin of attraction of the 
metastable CCW state decreases when decreasing \i. In 
fact, similar excited pulses and ISI distribution still hold 
for fi = 1.59 when the saddle and the metastable CCW 
state have disappeared. This can be understood by the 
fact that such dynamics arises as a scar of the bifurca- 
tions nearby |24j . 

All of the previous theoretical analysis has been done 
for the stochastic rate-equation system ([3])-([5]), where for 
reasons of simplicity the noise terms have been added 
in a symmetric way to both counter-propagating modes. 
From the experimental set-up explained in Section [TTJ 
one can wonder whether this symmetric assumption is 
valid as noise coming from the external SOA is injected 
mainly in one direction (the stable CW state). In order 
the check the validity of such an assumption, we have 
checked our analysis for asymmetric noise contributions, 
introducing noise with strength D cw = D in Eq. ([3| for 
the stable CW mode and noise with strength D ccw < D 
in Eq. Q for the metastable CCW mode. Qualitatively 
similar results are obtained as presented throughout this 



manuscript, which is reflected in Fig. 10 d) showing the 
ISI distribution for fi = 1.65, D cw = D = 2.5 x 1(T 4 
and different values of D ccw . 



-i 



V. DISCUSSION AND CONCLUSION 

In conclusion, we have investigated the features of 
excitability in optical systems which are close to Z2- 
symmetry and have elucidated the differences with other 



— more common — excitability scenarios [5 HTrJI |2"H 139] , 
The key message of this paper is the experimental and 
theoretical observation of a spread of the pulse amplitude 
and width in MRLs in the presence of noise, where a cor- 
relation between the amplitude and width of the excited 
pulses has been observed. Such a correlation between 
these two quantities is the signature of the deterministic 
evolution of the system once the separatrix is crossed. 
Multi-pulse excitability [24 is also present due to the 
finiteness of the noise intensity and the relaxation to the 
quiescent state being slow due to the closeness to a ho- 
moclinic bifurcation in parameter space. 

For a realization of an all-optical neural network based 
on MRLs, such a pulse width distribution can present a 
drawback compared to other excitable systems. However, 
such drawback is eventually balanced by the easy inte- 
grability of MRLs on an optical chip and the possibility 
to control the coupling between ring cavity and waveg- 
uide [43]. Furthermore, excitability in MRLs does not 
require feedback from an external cavity or optical injec- 
tion from a master laser as in the case of other optical 
systems. 
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